One-dimensional half-metallic interfaces of two-dimensional honeycomb insulators 
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We study zigzag interfaces between insulating compounds that are isostructural to graphene, 
specifically II- VI, III-V and IV-IV two-dimensional (2D) honeycomb insulators. We show that 
these one-dimensional interfaces are polar, with a net density of excess charge that can be simply 
determined by using the ideal (integer) formal valence charges, regardless of the predominant co- 
valent character of the bonding in these materials. We justify this finding on fundamental physical 
grounds, by analyzing the topology of the formal polarization lattice in the parent bulk materials. 
First principles calculations elucidate an electronic compensation mechanism not dissimilar to oxide 
interfaces, which is triggered by a Zener-like charge transfer between interfaces of opposite polarity. 
In particular, we predict the emergence of one dimensional electron and hole gases (IDEG), which 
in some cases are ferromagnetic half-metallic. 



Since the first two-dimensional material (graphene) 
was successfuly synthesized by mechanical exfoliation [T] , 
a growing number of studies has been devoted to other 
planar systems. These include BN monolayers [2], lay- 
ered transition metal oxides [5j , dichalcogenides [H |3] , or 
topological insulators such as Bi2Te3 or Bi2Se3 [5]. Of 
particular note in this context are a number of II- VI, 
III-V and IV-IV compounds, whose stable bulk phase 
is wurtzite, but may adopt a planar graphitic phase in 
ultrathin films. These were first predicted theoretically 
several years ago [71 [5] , and later grown in the lab [5] • 

Recent improvements in growth techniques, allow join- 
ing these nanosheets together. In addition to the many 
heterostructures studied so far involving graphene and 
BN, obtained by modified stacking [TU] or segregated 
nanosheets [TTl [Hj , there are and will be coplanar het- 
erostructures [13L I14j , made from different insulators on 
the same sheet many of them having a significant degree 
of ionic character. This brings up the important ques- 
tion of whether "polar discontinuities" [TS] might play a 
role in these low-dimensional systems. Vertical stacking 
of nanosheets is not dissimilar to oxide interfaces [H], 
and polarity effects will resemble those of ultrathin films 
already discussed in the literature [H |8]. However, to 
our knowledge, coplanar interfaces between 2D insulators 
have not been addressed from first principles, except for 
a qualitative discussion of the electrostatics [T7] . 

Compared with the case of polar oxide interfaces, there 
are a few additional challenges that are specific to the 2D 
case. First, the electrostatics is somewhat more compli- 
cated: unlike a plane of charge (interface between 3D 
materials) , which generates a uniform electric field in the 
whole space (essentially a ID problem), a line of charge 
(interface between 2D materials) produces a logarithmi- 



cally divergent potential with inhomogeneous stray fields 
in the vacuum region. Second, the symmetry group (with 
a main six-fold rotation axis but no center of symmetry) 
gives rise to different polarization classes than the ones 
found in centrosymmetric insulators [TH]. This calls for 
special care in the application of the modern theory of po- 
larization, and in particular of the interface theorem |19j , 
for the determination of the excess "bound" charge at a 
given interface. Moreover, the lack of center of symmetry 
implies that these materials are piezoelectrically active. 

Here we perform extensive density-functional theory 
calculations to show that (i) there indeed exists a net 
charge at the ID zigzag interface between two chemi- 
cally distinct polar insulating honeycomb domains; that 
(ii) the linear charge density is uniquely determined 
from bulk properties of the participating materials, and 
is related to the nontrivial topology of their respec- 
tive polarization lattice; that (iii) the divergent electro- 
static potential produced by the lines of charge triggers 
a Zener-like breakdown mechanism in the limit of in- 
creasingly wide domains; and that (iv) such mechanism 
gives rise to spin-polarized one- dimensional electron/hole 
gases (IDEG/IDHG). In what follows, we proceed ad- 
dressing points (i) to (iv) in the same sequential order. 

According to the modern theory of polarization |19) . 
and as further elucidated in Refs. [20l - f22| . the net charge 
at the interface between two compounds can be exactly 
given by invoking only bulk properties of the parent ma- 
terials. In particular, one needs to consider the lattice 
(Fig. [T]) of allowed values for the so-called formal polar- 
ization, P, of either crystalline constituent. Then, based 
on the interface orientation and termination j22j , one can 
readily deduce the interface net charge. An intuitive way 
to do this is by representing the electronic structure of the 
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material in terms of localized Wannier functions p51l24j . 
For polarization purposes, these effectively map the con- 
tinuous charge density distribution of the periodic crystal 
into a discrete set of classical point charges, located at 
the Wannier centers [SS] shown in Fig. [ija). Then the 
formal polarization of the crystal is simply given by the 
total dipole moment of the ensemble of these charges and 
the ion cores divided by the unit cell volume. Of course, 
there are infinite possible choices for the primitive basis 
of atoms and Wannier orbitals, hence the multivalued- 
ness of the formal P [see Fig. [T|b)]. 

A particularly convenient approach is that of com- 
bining the valence Wannier functions with their nearest 
neighbor ion core, obtaining a negatively charged object 
that we call "Wannier anion" [22^. This way, we ob- 
tain polarization lattices that include only a subset of the 
original points [larger symbols in Fig.[l|^b)]. More specif- 
ically, each of these restricted points simply corresponds 
to the dipole moment of a system of two integer point 
charges Q = ±iVe, placed at two arbitrarily chosen lat- 
tice sites, where N is the "formal valence" of the cation in 
each compound. Note that, due to the three-fold symme- 
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FIG. 1: (Color online) (a) The position of the Wannier centers 
(small circles) in the primitive unit cell of ZnO, AIN and SiC. 
The "cation" (Zn, Al, Si) is represented by squares, and the 
"anion" (O, N, C) by circles, "a" is the cell parameter (3.316 
A, 3.121 A, 3.110 A for ZnO, AIN and SiC respectively), and 
A is the distance of the sp'^ Wannier centers from the anion 
(A=0.431 A, 0.491 A, 0.677 A respectively), which increases 
with the bond covalency. (b) The triangular lattice of allowed 
polarization values for the II- VI (red circles), III-V (black tri- 
angles) and IV-IV (blue squares) 2D honeycomb compounds. 
Each symbol represents a polarization vector, indicated by 
P. The box highlights the values obtained with "Wannier 
anions" (see text), (c) Sketch of the zigzag interface between 
two compounds (AaBa in red and A^Bj, in blue). 



try of these structures, both the Wannier anion and the 
cation carry zero dipole moment; hence, the formal polar- 
ization only contains the point-charge contribution [22]. 
This has some analogies with the LaAlOa/SrTiOa case, 
where it was shown that covalency effects are also irrele- 
vant to charge-counting purposes [221 |5S] . 

Each arbitrary choice of atomic basis, in turn, yields 
a well-defined interface termination once the unit cell is 
periodically repeated to generate a semi-infinite domain. 
For the zigzag interfaces considered in this work, the ap- 
propriate choice of the bulk unit cell (i.e. one of the 
many choices that tile the whole system without leav- 
ing any unpaired ion at the interface) is that shown by 
shaded regions in Fig. [IJc). The dipole moment along y 
of such a unit cell is given hy d — —Qa/{2^/3), yielding 
a formal polarization of Py — —Q/3a. The interface the- 
orem relates differences in formal polarization to the net 
bound charge at the interface, and we have 
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This implicitly assumes that both parent compounds are 
hexagonal and have the same equilibrium lattice param- 
eter. In practice there is always a small lattice mismatch, 
and to realize a pseudomorphic interface one or both con- 
stituents need to be strained. Then, Eq. ([T]) must be 
corrected by piezoelectric contributions on Py. 

The electrostatic potential produced by this interfa- 
cial net charge diverges in the limit of increasingly wide 
domains and needs to be somehow neutralized. To inves- 
tigate possible compensation mechanism, we use explicit 
first-principles calculations of 2D superlattices, which are 
constructed by periodically repeating two composition- 
ally distinct nanoribbons along the stacking direction (y) , 
so that the interface direction (x) follows the zigzag edge 
of the hexagonal lattice [see Fig.jljc)], consistent with the 
above discussion of the bulk polarity. Having a relatively 
low (<7%) lattice mismatch, we chose SiC, AIN and ZnO 
as representative compounds for groups IV-IV, III-V and 
II- VI respectively. The periodicity of the superlattices is 
controlled by varying the nanoribbon width, m (given by 
the number of AB zigzag chains). Only stoichiometric 
compounds are considered, which means that two chem- 
ically distinct interfaces are present in our superlattices: 
one of type AaBaAbBb (A^: cation of material a, Bf,: an- 
ion of material b etc.) and one of type AbB^AaBa. Along 
the normal direction to the sheet we include a vacuum 
layer of 20A, sufficiently thick to obtain converged re- 
sults (no significant differences are observed with a vac- 
uum layer of up to 200A) . The calculations are performed 
using the spin-polarized Perdew-Burke-Ernzerhof (PBE) 
exchange correlation functional |26) as implemented in 
the Siesta code [27l |28| . Norm conserving pseudopoten- 
tials [ini and a double-^ polarized basis set are used to 
represent the valence electrons [301 • Atomic forces are 
relaxed to less than 20 me v/A. 
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FIG. 2: (Color online) (a) Spin-resolved layer-by-layer den- 
sity of states of the AlN-SiC (m=16) superlattice. The top 
half corresponds to the AIN slab, and the bottom half to the 
SiC slab, with the interface terminations indicated. Dashed 
(red) line schematizes the electrostatic potential, V(y), for 
an infinite array of alternating -I-/- cylindrical charges, (b) 
The spatial (plane averaged) charge distribution computed 
from the local density of states (LDOS) is shown in the lower 
panel (shaded curves are macroscopically averaged LDOS). 
(c) The density of free electrons (or holes) calculated as a 
function of slab width (m) for the two superlattice systems 
studied (AIN-SiC and ZnO-SiC). Theoretical values corrected 
by piezoelectric effects are marked by dashed lines. 



Figure [2]j a) shows the spin-resolved layer-by-layer den- 
sity of states (DOS) of the AIN-SiC (m=16) superlattice 
(III-V / IV- IV) . The presence of a macroscopic sawtooth- 
like potential, corresponding to electric fields of opposite 
polarity in the two material regions, is clearly visible from 
the plot. These fields are generated by the net interfacial 
charges discussed above. Under this potential the valence 
band maximum (VBM) at the p-type (Al-C) interface 



eventually becomes higher in energy than the conduction 
band minimum (CBM) at the n-type (Si-N) interface, 
and a charge transfer becomes favorable in energy. The 
resulting accumulation of free carriers along the inter- 
faces originate the aforementioned IDEG/IDHG whose 
spatial localization can be appreciated from their calcu- 
lated planar-averaged densities, plotted in Fig.[2]^b). The 
IDEG/lDHG, in turn, partially neutralize the net charge 
at either interface, reducing the potential offset. This 
feedback mechanism results in an effective pinning of the 
p-VBM and the n-CBM, which remain close in energy re- 
gardless of the value of m after breakdown occurs. Thus, 
many aspects of the present electronic mechanism resem- 
ble those already reported for oxide superlattices [25]. 
However, we stress that the electrostatics of this system 
is different, as evidenced by the macroscopic potential 
profile, V(y), followed by the DOS, which is analogous 
to that of an infinite array of alternating positive and 
negative linear charge densities [sketched in Fig.[2]ja)]. 

At smaller ribbon widths (not shown) the associated 
potential drop is smaller than the electronic band gap of 
either participating compound, and the system can re- 
main insulating (i.e. no "Zener breakdown" occurs). In 
the large m limit, on the other hand, the IDEG/lDHG 
densities tend to the "ideal" value, which corresponds 
to perfectly neutralizing the discontinuity in the nor- 
mal component of the bulk P, according to Eq. [T] This 
dependence of the charge transfer with ribbon width is 
shown in Fig. ^c) for both AIN-SiC (III-V / IV-IV) and 
ZnO-SiC (II- VI / IV-IV) interfaces. This nicely shows a 
cross-over scenario from a short-period uncompensated 
regime, to a longer-period regime where a Zener-like 
breakdown mechanism occurs, similar to the case of a 
LaAlOa-SrTiOs superlattice. The precise details of the 
cross-over depends on the band gaps and alignments be- 
tween the two materials, which might not be properly 
given by DFT calculations, although the physical mech- 
anism remains valid. Note that the asymptotic limits 
(dashed lines in Fig. [2jc)), in fact, do not correspond ex- 
actly to the ideal rational values (1/3 and 2/3) of Eq. [T] 
This is due to piezoelectric effects, which again are ac- 
counted for by macroscopic bulk properties of the par- 
ticipating materials. In particular, the lattice mismatch 
between SiC and ZnO, forces a compression in the ZnO 
ribbon (and a minor expansion in SiC) resulting in a shift 
of its formal polarization lattice along the axis normal to 
the interface, hence decreasing Psic — PznO ■ 

We find that the ID gases of carriers at the interfaces 
of these insulating materials are fully spin-polarized. To 
further investigate magnetic effects, we plot in Fig.[3]the 
band structure near the Fermi level, Ef, for a AIN-SiC 
superlattice with thickness that corresponds to the un- 
compensated regime (a) , and a partially compensated su- 
perlattice with larger thickness (panels (b)-(d) in the fig- 
ure). There are three dispersive tt bands close to Ef that 
are involved in the compensation mechanism, one occu- 
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pied and mostly localized at the Al-C interface, and two 
empty and mostly localized at the Si-N interface. The 
flat bands at X for the non-magnetic solutions Fig. [sj^a) 
and[3|^b) very much resemble the edge states in graphene 
nanoribbons described with a simple tight-binding model 
for the TT electrons [3T]. The corresponding high-density 
of electronic states at Ef suggests magnetic ordering due 
to a Stoner instability. 

Indeed, two spin-polarized phases with ferromagnetic 
orderings along the interface are lower in energy than the 
paramagnetic solution, and correspond to the transfer of 
charge (electrons) with well defined spin from one inter- 
face (p-VBM) to the other (n-CBM): when the spin is 
preserved in the transfer, as in Fig. [3][c), the electronic 
magnetization at both interfaces is antiferromagnetically 
aligned ( AF) , whereas a spin-flip of the transferred charge 
gives the ferromagnetic (FM) solution (panel d) . The two 
solutions, AF and FM between interfaces, are almost de- 
generate in energy, because the IDEG and the IDHG are 
decoupled by the wide insulating ribbon. Notice that the 
whole system is half-metallic for the AF solution, with 
full spin polarization of the conducting electrons. 

It is generally agreed that a one-dimensional metal is 
unstable with respect to a symmetry lowering modula- 
tion of the charge W'l . The coupling between lattice vi- 
brations and electrons near the Fermi level drives Peierls 
distortions that open a gap. Such deformation could also 
be a way to resolve the instability due to the large den- 
sity of states at the Fermi level shown in Fig. [sf^b) . To 
address the stability of the IDEG against Peierls dis- 
tortions, calculations using zone- folding techniques were 
performed for the AlN-SiC (m=16) superlattice, includ- 
ing up to three or four unit cells along the x axis (in the 
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limit of large m a supercell with three unit cells should 
be enough to accomodate symmetry-lowering deforma- 
tions capable of opening a gap in the right region of the 
Brillouin zone). We did not detect any trace of bond- 
length alternations in our tests, all of which remained in 
a metallic state. This is in agreement with previous stud- 
ies that showed that Peierls distortions can open a gap in 
similar edge bands only for narrow (<6A) graphene |33j 
and graphane nanoribbons [33]. 

In addition to competing structural/electronic ground 
states, another factor that could thwart experimental 
realization of the IDEG/IDHG predicted here are de- 
fects. In particular, these polar structures may, in cer- 
tain growth conditions, find more energetically favorable 
compensation mechanisms than electronic transfer. For 
example an anionic vacancy could donate electrons to 
the n-type interface, or combinations of anionic vacancies 
and cationic vacancies at each interface could provide the 
necessary charge compensation without free electrons. 
Similar mechanisms have been proposed at the LaAlOa- 
SrTiOs interface and in ferroelectric thin films [35 l 136 ) . 
Substitutional defects and charge transfers between edge 
states in isostructural C/BN superlattices have been also 
discussed 37J. If the growth conditions can be controlled 
such that these point defects arc minimized it may be 
possible to utilize the systems here for various applica- 
tions such as spintronics, sensors and electronics. 

In conclusion, we have shown that polar compensa- 
tion mechanisms can give rise to fully spin-polarized one- 
dimensional electron and hole gases at the interface be- 
tween two coplanar insulating bidimensional materials, in 
close analogy with the 3D case of LaAlO3/SrTiO3(001) 
interface. The net charge at the interface can be deter- 
mined solely from bulk properties (the polarization) of 
the parents materials, as required by the interface theo- 
rem [19] . The polarization is given by a simple charge 
counting procedure, irrespective of covalency/ionicity, 
that results in three distint sets of allowed polarization 
lattices, for the II- VI, III-V and IV-IV compounds, and 
to the formation of new electronic states (IDEGs) at the 
interface between two materials of different classes. 
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FIG. 3: (Color online) The band structure of the AIN-SiC 
superlattices for m = 8 (a), and m = 16 with paramagnetic- 
NM (b), antiferromagnetic-AF (c) and ferromagnetic-FM (d) 
solutions. Solid (black) and dashed (red) lines correspond to 
bands with different spin components. The inset in (b) shows 
the Fermi level crossing of bands at the zone boundary. 
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